Strong renormalization of the Fermi-surface topology close to the Mott transition 
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The underlying Fermi surface is a key concept for strongly-interacting electron models and has 
been introduced to generalize the usual notion of the Fermi surface to generic (superconducting or 
insulating) systems. By using improved correlated wave functions that contain backfiow and Jastrow 
terms, we examine the two-dimensional t—t' Hubbard model and find a non-trivial renormalization 
of the topology of the underlying Fermi surface close to the Mott insulator. Moreover, we observe a 
sharp crossover region, which arises from the metal-insulator transition, from a weakly interacting 
metal at small coupling to a resonating valence-bond superconductor at intermediate coupling. A 
violation of the Luttinger theorem is detected at low hole dopings. 

PACS numbers: 71.27.+a,71.18.+y,71.30.+h,74.20.-z,74.72.-h 



I. INTRODUCTION 

The single-band Hubbard model with extended hop- 
ping on the square lattice has been widely investigated 
since the appearance of the high-temperature supercon- 
ductors. Indeed, it is believed to represent the minimal 
model that is necessary to describe the electronic corre- 
lations in the Copper-Oxygen planes of Cuprate materi- 
als. A very rich variety of phases has been discussed in 
this context, including antiferromagnetism, superconduc- 
tivity, charge-density waves, and non-Fermi-liquid met- 
als i Studies with various numerical techniques, rang- 
ing from dynamical mean- field theory (DMFT) includ- 
ing its cluster extensions^— to quantum Monte Carlo 
techniques, 8 — as well as analytic approaches ) 13 ! 14 have 
addressed hotly debated topics like the nature of the 
pseudo-gap phase or the superconducting correlations. 

Landau 1 ^ and Luttinger— have shown that the Fermi 
surface, which is the locus of gapless electronic excita- 
tions in fc-space, represents a pivoting concept in the 
theory of Fermi liquids. The generalization of this no- 
tion for gapped systems (e.g., superconductors or Mott 
insulators) leads to the idea that there is an underlying 
Fermi surface that becomes gapped because of some sym- 
metry breaking (leading to superconductivity) or elec- 
tronic correlation (leading to a Mott insulator). In par- 
ticular, the underlying Fermi surface can be defined by 
the locus of points where ReG(k,uj = 0) changes sign 
(G(k,uj) is the single-particle Green's function), pass- 
ing either to infinity (for usual Fermi liquids with well 
defined quasiparticles) or to zerOfi 7 - The question of de- 
termining the underlying Fermi surface is of central im- 
portance in strongly correlated systems, particularly in 
view of experiments with angle-resolved photoemission 
spectroscopy (ARPES).— From the theoretical point of 
view, there are few attempts to study the topology of the 
Fermi surface in correlated systems; furthermore, they 
are limited to mean-field approaches,— ~— including re- 
cent calculations with cluster DMFT^ Here, we make 



a substantial step forward and consider non-perturbative 
calculations directly in a two-dimensional system. We 
use variational wave functions containing backfiow cor- 
relations, which have been shown to be very accurate in 
the Hubbard model ) 23 ' 24 and focus our attention on the 
topology of the underlying Fermi surface and the super- 
conducting properties. We show that a strong renormal- 
ization of the underlying Fermi surface takes place close 
to half filling, when the Mott insulator is approached. In 
particular, the renormalization to perfect nesting occur- 
ring at the interaction-driven metal-insulator transition 
opens a new perspective on the crucial role of momentum 
dependence in describing the Mott transition. 

At finite dopings we find a sharp crossover line, which 
emerges from the Mott-Hubbard transition point at 
half filling and separates a weakly-correlated metal and 
a strong-coupling superconducting state, in agreement 
with a recent observation of an unconventional metal- 
lic state at finite dopings^ Our findings could be also 
related to a recently proposed first-order line separat- 
ing two metallic phases, one with a pseudogap and one 
without^ For the superconducting order parameter we 
obtain a sizable signal for moderately large on-site inter- 
actions, i.e., U/t ~ 7— 10. 

The paper is organized as follows: in section UH we 
introduce the Hamiltonian and describe our variational 
wave function; in section Mil we present our numerical 
results and, finally, in section lTVl we draw the conclusions. 



II. MODEL AND VARIATIONAL WAVE 
FUNCTION 

We consider the Hubbard model with extended hop- 
ping on a two-dimensional square lattice, 

H = -tJ24 a c ja -t' ]T4c iCT +h.c.+t/^n it n u , (1) 

(ij)a- » 
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FIG. 2: (Color online) Variational hopping ratio t'/t as a 
function of doping 8 and U/t for the case with t'/t — —0.4. 
The density plot is obtained by using results with L = 162 
and 242. 



FIG. 1: (Color online) Upper panel: variational hopping ratio 
t'/t as a function of U/t for t'/t = -0.4 and -0.75. Full 
(empty) symbols refer to the presence (absence) of backflow 
correlations in the variational wave function. Lower panel: 
variational hopping ratio t'/t as a function of the doping 5 
for U/t = 5 and 10 with t'/t — —0.4. Here, only results with 
backflow terms are shown. Data are shown for a L = 162 
lattice; in the lower panel few points on a L = 242 lattice are 
added close to half filling. 

where c\ a (c ia ) denotes the electron creation (destruc- 
tion) operator of one electron on site i with spin a = J [, J., 
(ij) and ((ij)) indicate nearest and next-nearest neighbor- 
sites respectively; rii a — ^ ia c ia is the electron density; t 
and t' are the nearest and next-nearest neighbor hopping 
amplitudes, and U is the on-site Coulomb repulsion. Cal- 
culations are performed on 45-degree tilted clusters with 
L = 2 x I 2 sites (I being an odd integer) and periodic 
boundary conditions. The number of electrons is N, such 
that the hole doping is 5 = 1 — n, with n = N/L. 

The question of determining the underlying Fermi sur- 
face has been addressed by using a renormalized mean- 
field approach in Refs. [l9ll20l Here, we include elec- 
tronic correlations in a non-perturbative way. In a first 
step, we construct uncorrelated wave functions given by 
the ground state |BCS) of a superconducting Bardeen- 
Cooper-Schrieffer (BCS) Hamiltonian^^ 

Hbcs = £ betray + £ A fe 4 t cL H + h.c, (2) 

ka k 

where both the free-band dispersion and the pairing 
amplitudes are variational functions. We use the 
parametrization 

= — 2t(cos k x + cos k y ) — At' cosk x cosk y — fi (3) 
A k = 2A B cs(cosfc i: - cosfcy), (4) 



where the effective hopping amplitude t' , the effective 
chemical potential /i, and the local pairing field Abcs are 
variational parameters to be optimized. The parameter 
i is kept fixed to set the energy scale. We also consid- 
ered longer-range effective hopping parameters in Eq. (U) 
finding that all hoppings beyond the ones present in the 
Hamiltonian ([TJ are optimized to zero. 

The correlated state I^bcs)^ without backflow terms, 
is then given by |*bcs) = J|BCS), where J = 
exp(— 1/2 y^-j VijUiUj) is a density-density Jastrow factor 
(including the on-site Gutzwillcr term vu), with the uy's 
being optimized for every independent distance \i — j\. 
Notably, within this kind of wave function, it is possi- 
ble to obtain a pure (i.e., non-magnetic) Mott insulator 
for a sufficiently singular Jastrow factor v q ~ 1/q 2 (y q 
being the Fourier transform of Vy), while a supercon- 
ducting (metallic) state is found whenever v q ~ 1/q and 
Abcs > (A BC s = 0)ja 

A size-consistent and efficient way to further improve 
the correlated state I^bcs) for large on-site interactions 
is based on backflow correlations. In this approach, each 
orbital that defines the unprojected state |BCS) is taken 
to depend upon the many-body configuration, such to in- 
corporate virtual hopping processes. 2 ^ 24 All results pre- 
sented here are obtained by fully incorporating the back- 
flow corrections and optimizing individually every varia- 
tional parameter in and A/., in the Jastrow factor J7\ 
as well as for the backflow corrections. 



III. RESULTS 

A. Fermi-surface renormalization 

We start our analysis by discussing the evolution of the 
underlying Fermi surface topology as a function of hole 
doping 5 and interaction strength U/t. Within the varia- 
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FIG. 3: (Color online) Static spin-spin correlations S(q) as a 
function of \q\/ir along the F = (0, 0) to M — (n, ir) direction, 
for t'/t = —0.4. Data are presented for U/t = 5 (squares), 
U/t — 5.6 (triangles), U/t = 6 (circles) and U/t = 7 (dia- 
monds). Full (empty) symbols refer to L = 242 (L = 162). 
The location of the point Q = (ir, it) is marked with an arrow. 



tional approach, the underlying Fermi surface can be eas- 
ily defined and corresponds to the highest occupied mo- 
menta Indeed, Ef- = v£fc~+ ^'t corresponds, within 
renormalized mean-field theory^ to the excitation spec- 
trum of projected Bogoliubov modes and hence to the 
dispersion of the renormalized quasiparticles. Here, we 
would like to stress that, although corresponds to the 
"non-interacting" spectrum, it may be strongly renor- 
malized because of the full optimization of the variational 
wave function: due to the presence of Jastrow and back- 
flow terms, non-trivial values of t and Abcs may be 
obtained. 

Firstly, we consider the half-filled case, where, by in- 
creasing the ratio U/t a metal-insulator transition is en- 
countered. We present the results for t'/t = —0.4 and 
t'/t = -0.75, see Fig. HJ In both rcnormaliza- 
tion towards perfect nesting, namely t'/t = 0, is shown 
to occur both at the metal-insulator transition and in 
the limit U/t — » oo. In the metallic phase, the ratio t'/t 
is only slightly modified with respect to its bare value, 
and this result is found both with and without back- 
flow correlations, demonstrating that even the simple 
wave function without backflow may capture a correct 
description of the metallic phase. Then, at the metal- 
insulator transition, the variational ratio of the hopping 
parameters is strongly renormalized in presence of back- 
flow correlations, i.e., t'/t — >• 0, driving the underlying 
Fermi surface to be perfectly nested. This fact suggests 
that a full momentum resolution is a crucial ingredient 
to properly describe a metal-insulator transition. Re- 
markably, this renormalization does not occur if back- 
flow correlations are not included. By further increasing 
the on-site Coulomb repulsion, the renormalized hopping 
ratio changes sign and then decreases again to zero at 
large U /t£H This renormalization, which is again possi- 
ble only in presence of backflow correlations, is in agree- 
ment with renormalized mean- field studies of the t—J 
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FIG. 4: (Color online) Left panels: the underlying Fermi 
surface defined by £k = (solid blue lines) and the non- 
interacting Fermi surface £° = (dashed black lines) for the 
same parameter values of the corresponding right panels (see 
text for more details). Right panels: momentum distribution 
function rtk along the path in the Brillouin zone connecting 
the points F = (0,0), M = (ir,ir), and X = (n,0). Data are 
shown at U/t = 10 and t'/t = -0.4 for a L = 242 lattice size 
at three dopings (S — 0, 0.033, and 0.174). Arrows indicate 
the position of the underlying Fermi surface, as found in left 
panels. 



model, 20 where the Fermi surface renormalizes to per- 
fect nesting for 5 — > 0, and with a slave-spin study of 
the Hubbard model in the large-interaction limit<21 Note 
that for U/t > 25 the numerical accuracy for the ratio 
t'/t is not enough to distinguish a finite from a vanishing 
value, see Fig. Q] 

We would like to mention that a Fermi-surface renor- 
malization close to the Mott transition has already been 
described in the one-dimensional Hubbard models Also 
in that case, a perfectly nested Fermi surface, i.e., kp = 
±7r/2, is found. This consideration enforces the idea that 
the variational hopping ratio in two dimensions does not 
vanish accidentally at the Mott transition, but is renor- 
malized in order to get perfect nesting of the underlying 
Fermi surface. 

The variational hopping changes very rapidly by dop- 
ing the Mott insulator and tends towards the bare value 
at high values of hole dopings, see the case U/t = 10 
in Fig. [TJ In summary, we present in Fig. [5] a complete 
density plot for the variational hopping t'/t as a func- 
tion of doping 5 and U/t, with bare t'/t = —0.4. While 
we have almost no renormalization when the half-filled 
case is metallic (i.e., for U/t < 6), we observe that a 
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FIG. 5: (Color online) Left panel: Ratio between the number 
of electrons enclosed by the renormalized Fermi surface, h, 
and the number of electrons effectively present in the system, 
n, as a function of doping. Data refer to the case U/t = 10 
and L = 242. Right panel: Same quantity as in the left panel 
as a function of U/t at half filling, for L — 242. 



strong renormalization close to half filling is present for 
U/t>6. Remarkably, the metal-insulator transition oc- 
curring at half filling comes with a sharp crossover line at 
finite doping that we can characterize through an abrupt 
change in the value of i'/i. This crossover separates a 
weakly-correlated metal at low U/t and a strongly renor- 
malized state at intermediate/strong values of U ft. This 
result is in agreement with a recent ARPES study on 
YBa2Cu306+ XJ which found indications for an uncon- 
ventional metallic state. 25 



B. Magnetic correlations 

We provide here further evidence of the sharp crossover 
line by calculating the static spin-spin correlations, de- 
fined as 



-Rn) 



(5) 



where is the z-component of the spin operator on 
site to. The presence of (short-range) antiferromagnetic 
correlations is signaled by the appearance of a (non- 
diverging) peak in S(q), located at Q = (tt,tt). As 
shown in Fig. [3] for the doping S = 0.025, the corre- 
lated resonating- valence bond (RVB) state at U/t > 6 is 
characterized by antiferromagnetic correlations that are 
strongly enhanced with respect to the metallic phase at 
U/t < 6. The two regimes are clearly separated by a 
jump in the value of S(Q); the short-range nature of the 
antiferromagnetic correlations is confirmed by a size scal- 
ing study. 



C. Momentum distribution function 

In Fig. [2J we compare the results of the underlying 
Fermi surface with the ones for the momentum distri- 
bution nk = (c].<T C ka)' f° r = 10 and three differ- 
ent dopings. Interestingly, in the low-doping regime, 
the area enclosed by = is different from the one 
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FIG. 6: (Color online) Pairing amplitude AbcsA as a func- 
tion of doping for t' /t = —0.4 and different values of U/t. 
Data are presented for a L = 162 lattice size. 



enclosed by the non-interacting system £° = 0, where 
= — 2i(cos k x + cos k y ) — At' cos k x cos k y — fi , Mo being 
the bare chemical potential. At half filling, the underly- 
ing Fermi surface is closed, but the system is insulating 
and consistently shows a completely smooth behav- 
ior. This is due to the presence of a strong Jastrow fac- 
tor v q ~ 1/q 2 which is able to remove the singularities 
present in |BCS). 32 As soon as a small doping is consid- 
ered, the system becomes conducting and a finite jump 
in rife appears along the nodal direction T —¥ M , where 
the pairing amplitude A^ vanishes. At small dopings, the 
variational hopping t ft undergoes an abrupt change (see 
Fig. [2]) and, therefore, the Fermi surface becomes rapidly 
open, see Fig. S) Finally, at large dopings, the supercon- 
ducting gap vanishes and a full Fermi surface is recovered 
also in the momentum distribution, i.e., a finite jump is 
also detected along the M — > X direction. Additional 
results on the evolution of the underlying Fermi surface 
as a function of doping for U/t = 10 and U/t — 5 can be 
found in the Appendix [A"l 



In the large doping case, the jump along the M — > X 
direction coincides with the position of the underlying 
Fermi surface. When the doping is small and there is 
no jump along the M — s> X direction, one could extend 
the previous concept and associate the location of the 
maximal gradient of with the position of the under- 
lying Fermi surface. However, as shown in Fig. 01 this 
approach would lead to incorrect results, since the under- 
lying Fermi surface is closer to the X point with respect 
to the location of the maximal gradient. Therefore, our 
results suggest that, when doping is small and electronic 
correlation is important, some caution should be taken 
in deriving the Luttinger surface from the location of the 
maximal gradient of nfc ; 33 i 34 similarly to what has been 
discussed in Ref. [20 for the t— J model. 
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FIG. 7: (Color online) Superconducting order parameter squared 2 as a function of doping for t'/t = —0.4 and U/t = 7 (left 
panel) and U/t = 10 (right panel). 



D. Luttinger sum rule 

The renormalization of the underlying Fermi surface 
leads to a violation of the Luttinger sum rule. Let us 
denote by n the number of electrons that are enclosed 
by the renormalized Fermi surface = 0. As shown in 
Fig. SI close to half filling, the = contour is electron- 
like, containing less electrons than n, i.e., n/n < 1, while 
for 0.05 < 8 < 0.2, the contour is hole-like and n/n > 1. 
For larger values of the doping, the Luttinger count holds. 
The summary for the Luttinger sum rule is reported in 
Fig. [5] Our results are slightly different from the ones 
found at the simple mean-field level^ especially close to 
the Mott insulator where we find that Sn ~ — 0.1. In 
Fig. [SJ we also present the violation of the Luttinger sum 
rule at half filling, as a function of U/t. In the insulating 
state (U/t > 6) the underlying Fermi surface defined by 
£fc = encloses less electrons than n. This result is com- 
patible with the electron-like nature of the contour, due 
to the strong renormalization of the variational hopping 
parameters. On the contrary, in the metallic region we 
only observe a tiny violation, in particular close to the 
metal-insulator transition. This tiny violation can be re- 
lated to the amount of correlation present in the \ow-U 
metallic state and it is indeed significantly smaller than 
the degree of correlation that can be observed in the RVB 
state emerging from the Mott insulator upon doping, see 
FigEJ 



E. Pairing amplitude and pair-pair correlations 

We discuss now the variational pairing amplitude 
Abcs , which measures the tendency to create resonating 
singlets, and the actual pair-pair correlations (A(r)) = 
(S r So), where S$ = c^cj +a .^ - c^<4 +2 , t . At weak- 
coupling (when a conducting state is found at half fill- 
ing) , Abcs remains relatively small and does not change 
much with doping, see Fig [6] For U/t > 6 (when an in- 
sulating state is found at half filling) , Abcs has a sizable 
value and decreases monotonically as a function of the 
doping. This behavior of the pairing amplitude has been 



related to the pseudo-gap of the normal phase of Cuprate 
materials^ 

In Fig [JJ we report the results for the superconduct- 
ing order parameter <f> 2 = limj._j.oo A (r), similarly to 
what has been done in previous calculations for the Hub- 
bard and t — J modelsi 10 ' 35 Our results are presented 
for U/t = 7 and 10, where superconductivity develops at 
finite doping. For smaller values of the on-site interac- 
tions, a much smaller signal is obtained. By increasing 
U/t, the optimal doping becomes larger and the magni- 
tude of the order parameter decreases, suggesting that 
intermediate values of U/t are optimal for maximal crit- 
ical temperatures T c , when assuming that T c scales with 



IV. CONCLUSIONS 

Our calculations represent a first attempt to trace the 
Fermi surface renormalization in a truly two-dimensional 
system where the electron-electron correlation is treated 
beyond simple mean-field approaches. A substantial de- 
viation from the Luttinger sum rule is observed for large 
values of U/t, where the topology of the underlying Fermi 
surface changes from electron-like to hole-like by increas- 
ing the doping. It would be very interesting to verify 
these results on correlated materials, like Cuprate su- 
perconductors, where ARPES probes should be able to 
detect changes in the topology of the Fermi surface. 

Furthermore, we show evidence of a sharp crossover re- 
gion that originates from the metal-insulator transition 
and separates a weakly-correlated metal from a more cor- 
related RVB superconductor at low dopings. 

L.F.T. and C.G. acknowledge the support of the Ger- 
man Science Foundation through the Transregio 49. 



Appendix A: Underlying Fermi surface 

In Fig. [51 we present a systematic plot of the underlying 
Fermi surface defined by the contour = and of the 
non-interacting Fermi surface, £° = 0, as a function of 



6 



doping for the case U/t — 10. The underlying Fermi 
surface evolves from electron-like, when 8 < 0.05, to hole- 
like, for 5 > 0.05. When 8 ~ 0.2 the underlying Fermi 
surface becomes very close to the non-interacting one, 
as expected when doping becomes large and electronic 
correlation consequently less important. 

In Fig. [SJ we present also a plot of the underlying Fermi 



surface as a function of doping for U/t = 5. In this 
regime, the half- filled case is metallic and the renormal- 
ization of the variational parameters is very weak. Con- 
sequently, the underlying Fermi surface does not show a 
remarkable evolution as a function of doping and almost 
coincides with the non- interacting contour. 
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FIG. 8: (Color online) The underlying Fermi surface defined by = (solid blue lines) and the non-interacting Fermi surface 
= (dashed black lines) as a function of doping. Data are presented at U/t = 10 for L — 242. 




FIG. 9: (Color online) The underlying Fermi surface defined by = (solid blue lines) and the non-interacting Fermi surface 
£2 = (dashed black lines) as a function of doping. Data are presented at U/t — 5 for L — 242. 



